Rashba spin splitting in two-dimensional electron and hole systems 
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In two-dimensional (2D) hole systems the inversion asymmetry induced spin splitting differs 
remarkably from its familiar counterpart in the conduction band. While the so-called Rashba spin 
splitting of electron states increases linearly with in-plane wave vector k« the spin splitting of heavy 
hole states can be of third order in fcy so that spin splitting becomes negligible in the limit of small 
2D hole densities. We discuss consequences of this behavior in the context of recent arguments on 
the origin of the metal-insulator transition observed in 2D systems. 
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At zero magnetic field B spin splitting in quasi two- 
dimensional (2D) semiconductor quantum wells (QW's) 
can be a consequence of the bulk inversion asymmetry 
(BIA) of the underlying crystal (e.g. a zinc blende struc- 
ture) and of the structure inversion asymmetry (SIA) r 
the confinement potential. This B = spin splitting 
is the subject of considerable interest because it concerns 
details of energy band structure that are important in 
both fundamental research and electronic device applica- 
tions (Refs. p|-pl| and references therein). 

Here we want to focus on the SIA spin splitting which 
is usually, the dominant part of B = spin splitting in 2D 
systemsE To lowest order in ku SIA spin splitting in 2D 
electron systems is given by the so-called Rashba model,cl 
which predicts a spin splitting linear in ku. For small 
in-plane wave vector k\\ this is in good agreement with 
more accurate numerical computations .□ For 2D hole sys- 
tems, on the other hand, the situation is more compli- 
cated because of the fourfold degeneracy of the topmost 
valence band Tg, and so far only numerical computa- 
tions on hole spin splitting have been performed.trD In the 
present paper we will develop an analytical model for the 
SIA spin splitting of 2D hole systems. We will show that 
in contrast to the familiar Rashba model the spin split- 
ting of heavy hole (HH) states is basically proportional 
to k?,. This result was already implicitly contained in 

several numerical computations.ifl But a clear analytical 
framework was missing. We will discuss consequences of 
this behavior in the context of recent arguments on the 
origin oLthe metal-insulator transition observed in 2D 
systems E l 

First we want to review the major properties of the 
Rashba modelcl 



fff c ° = akxE-<7. 

In this equation a — (a x , a y , a z ) denotes thi 
matrices, a is a material-specific prefactor,F ' 1 



(1) 



uli spin 
and E 



is an effective electric field that results from the built-in 
or external potential V as well as from the position de- 
pendent valence band edge. For E = (0, 0, E z ) Eq. (|l|) 
becomes (using explicit matrix notation) 
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by analyzing 



obtain for the spin splitting of the energy dispersion 

f 6 s c °(k||) = ±(«^)fc|| (3) 

where k|| = (k x ,k y ,0). Using this simple formula sev- 
eral groups determined the prefactor 
Shubnikov-de Haas (SdH) oscillations .EHf 

Equation (^) predicts an SIA spin splitting which is 
linear in fc|| . For small k\\ Eq. (|^) thus becomes the dom- 
inant term in the energy dispersion £±(k||), i.e., SIA spin 
splitting of electron states is most important for small 2D 
densities. In particular, we get a divergent van Hove sin- 
gularity of thp density-of-states (DOS) at the bottom of 
the subbandla which is characteristic for a k linear spin 
splitting. As an example, we show in Fig. [l] the self- 
consistently calculated^ subband dispersion £±(k\\), DOS 
effective mass m* /mo, and spin splitting £+(kn ) — £-(k\\ ) 
for an MOS inversion layer on InSb. For small fcy the 
spin splitting increases linearly as a function of ku , in 
agreement with Eq. (JsJ) . Due to nonparabolicity the,spin 
splitting for larger ku converges toward a constants 

The spin splitting results in unequal populations N± 
of the two branches £±(kn). For a given total density 
N s = N + + iV_ and a subband dispersion £±(fcii) = 
± (aE z )ku with \i = h 2 /2m* we obtain 



N+ = -N„ ± 



(aE z ) 
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with k± = k x ±ik y . By means of perturbation theory we 



This equation can be directly pjcpmpared with, e.g. 
results of SdH experiments. E3~Ej 

The Rashba model (jl]) can be derived by purely group- 
theoretical means. The electron states in the lowest 
conduction band are s like (orbital angular momentum 
1 = 0). With spin-orbit (SO) interaction we have total 
angular momentum j = 1/2. Both k and E are polar 
vectors and k x E is an axial vector (transforming ac-. 
cording to the irreducible representation of T^O'Ej 
Likewise, the spin matrices <j x .<j v , and a z form an axial 
vector cr. The dot product (Q) of k x E and cr there- 
fore transforms according to the identity representation 
Ti, in accordance with the theory of invariants of Bir 
and PikusBj In the Tg conduction band the scalar triple 
product (Q) is the only term of first order in k and E that 
is compatible with the symmetry of the band. 
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Now we want to compare the Rashba model (|l|) with 
the SIA spin splitting of hole states. The topmost va- 
lence band is p like (I = 1). With SO interaction we 
have j = 3/2 for the HH/LH states (T v s ) and j = 1/2 for 
the SO states (T7). For the Tg valence band there are 
two sets of matrices which transform like an axial vec- 
tor namely J = (J x , J y , J z ) and J = ( J%, J^, Jf ) (Refs. 
^2| , p4] ). Here J^, J y and </ z are the. angular momentum 
matrices for j = 3/2. Thus we getcJ 

flf° = ft k x E • J + fa k x E • J. (5) 

Similar to the Rashba model the first term has axial sym- 
metry with the symmetry axis being the direction of the 



electric field E. The second term is anisotropic, i.e., it 
depends on both the crystallographic orientation of E 
and k. Using k • p theory we find that the prefactor /?2 is 
always much smaller than f3\, i.e., the dominant term in 
Eq. (||) is the first term. This can be easily understood 
by noting that the k • p coupling between Tg and Tg is 
isotropic, so that it contributes to /3i but not to /?2- The 
prefactor 02 stems from k • p coupling to more remote 
bands such as the p antibonding conduction bands T| 
and T c 7 . 

For E = (0, 0, E z ) Eq. (j|) becomes (using explicit 
matrix notation with j = 3/2 eigenstates in the order 
U = +3/2, +1/2, -1/2, -3/2) 
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fV3fe + 5/2fe_ 

5/2fc+ f%/3fe_ 

\ 3/4 fc_ lV3k+ / 



(6) 



Here the first term couples the two LH states (j z = 
±1/2), and it couples the HH states (j z — ±3/2) with 
the LH states. But there is no fe linear splitting of the 
HH states proportional to (3\ . The second matrix in Eq. 
(||) contains a k linear coupling of the HH states. 

We want to emphasize that iff,? and are effec- 

tive Hamiltonians for the spin splitting of electron and 
hole subbands, which are implicitly contained inrthe. full 
multiband Hamiltonian for the subband problemBO 



H = ff k . p (k||,fc 2 = -id z ) + eE z zt. 



(7) 



Here i/k-p is a k • p Hamiltonian for the bulk band 
structure (i.e., ffk p does not contain or Hjj®) 

and we have restricted ourselves to the lowest order 
term in a Taylor expansion of the confining potential 
V(z) = V + eE z z + 0(z 2 ) which reflects the inversion 
asymmetry of V(z). The effective Hamiltonians W) and 
(^) stem from the combined effect of -f/k-p and the term 
eE z z. For a systematic investigation of the importance 
of the different terms in H we have developed a novel, 
analytical approach based on a perturbative diagonaliza- 
tion of H using a suitable set of trial functions and using 
Lowdin partitioning.E3c3 Though we cannot expect ac- 
curate numerical results from such an approach it is an 
instructive complement for numerical methods, as we can 
clearly identify in the subband dispersion £(kii ) the terms 
proportional to E z which are breaking the spin degener- 
acy. Neglecting in i?k p remote bands like Tg and Ty we 
obtain for the SIA spin splitting of the HH states 



4°(fc||)cx±(/3 1 £ z )fef. 



(8a) 



In particular, we have no k linear splitting (and P2 = 
if we restrict ourselves to the Luttinger HamiltonianE-3 
which includes Tg and Ty by means of second orde* 
perturbation theory.El Accurate numerical computationsu 
show that the dominant part of the k linear splitting of 
the HH states is due to BIA. However, for typical den- 
sities this k linear splitting is rather small. For the LH 
states we have 



££g(fc||) « ±(PiE,)kn 



(8b) 



Thus we have a qualitative difference between the spin 
splitting of electron and LH states which is proportional 
to feii and the splitting of HH states which essentially is 
proportional to fc|. The former is most important in the 
low-density regime whereas the latter becomes negligible 
for small densities. Note that for 2D hole systems the 
first subband is HH like so that for low densities the SIA 
spin splitting is given by Eq. ( pa] ) . In Eq. (JsJ) the lengthy 
prefactors depend on the details of the geometry of the 
QW. Moreover, we have omitted a weak dependence on 
the direction of k|| . But the order of the terms with re- 
spect to feii is independent of these details. It is crucial 
that, basically, we have 
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with Ao the SO gap between the bulk valence bands 
and r", i.e., we have no SIA spin splitting for Ao = O.Ei 
This can be most easily seen if we express -ffk-p in a basis 
of orbital angular momentum eigenstates. 

A more detailed analysis of our analytical model shows 
that both Hq® and flf ° stem from a third order pertur- 
bation theory for k±, k z — —id z , and eE z z. This seems 
to be a rather high qeder. Nevertheless, the resulting 
terms are fairly large. Ell In agreement with Refs. p|J^Jl7| 
this is a simple argument to resolve the old controversy 
based on an argument by Andon3 that spin splitting in 2D 
systems ought to be negligibly small because for bound 
states in first order we have (E z ) = 0. We note that 
the present ansatz for the prefactors a and fli , 02 is quite 
different from the ansatz in Ref. [lO] We obtain Hq® and 
Ha® by means of Lowdin partitioning of the Hamiltonian 
"]) whereas in Ref. [lO] the authors explicitly introduced 
into their model. Moreover, we evaluate the matrix 
elements of eE z z with respect to envelope functions for 
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the bound states whereas in Rcf. 10 the authors consid- 



ered matrix elements of eE z z with respect to bulk Bloch 
functions. The latter quantities are problematic because 
they depend on the origin of the coordinate frame. 

As an example, we show in Fig. |^ the self-consistently 
calculatedO anisotropic subband dispersion £±(k||), DOS 
effective mass m* /mo, and spin splitting f+(kii) — £_(kii) 
for a [001] grown GaAs/ALj.sGao.sAs heterostructure. 
The calculation was based on a 14 x 14 Hamiltonian (r§, 
rf , rg, rg, and r?). It fully took into account both SIA 
and BIA. The weakly divergent van Hove singularity of 
the DOS effective mass at the subband edge indicates 
that the k linear splitting is rather small. (Its dominant 
part is due to BIA.E3) Basically, the spin splitting in Fig. 
H is proportional to fcjj. 

Only for the crystallographic growth directions [001] 
and [111] the hole subband states at ku = are pure HH 
and LH states. For low-symmetry growth directions like 
[113] and [110] we have mixed HH-LH eigenstates even 
at fc|| — 0, though often the eigenstates can be labeled by 
their dominant spinor components O The HH-LH mixing 
adds a k linear term to the splitting ( |8a| ) of the HH states, 
which often exceeds (/32E z )k\\ . However, this effect is still 
small when compared with the cubic splitting. 

For a HH subband dispersion £±(k^) = (fj)ku ± 

(PiE z )k^ we obtain for the densities N± in the spin-split 
subbands 



N ± = iiV s ± 



with 



X = 1 



1 - 4tt N s 



(PiE z 



(10b) 



The spin splitting according to Eq. ( |l0| ) is substantially 
different from Eq. (|J). For electrons and a fixed electric 
field E z but varying N s the difference AA^ = N + — 2V_ 

1 /2 

increases like N s whereas for HH subbands it increases 

3 /2 

like N s ' . Using a fixed density N s but varying E z it is 
more difficult to detect the difference between Eqs. (^) 
and ([To]). In both cases a power expansion of AN gives 
A A" = ai \E z \ +a 3 |£; z | 3 + C(|£; Z | 5 ) with a 3 < for elec- 
trons and a 3 > for HH subbands. 

The proportionality is completely analogous to the. 
effective g factor in bulk semiconductors EJ LassnigEZI 
pointed out that the B = spin splitting of electrons can 
be expressed in terms of a position dependent effective 
g factor g*(z). In the following we want to discuss the 
close relationship between Zeeman splitting and B = 
spin splitting from a more general point of view. Note 
that in the presence of an external magnetic field B we 
have k x k = (— ie/7i)B and the Zeeman-splitting in the 



rg conduction band can be expressed aa 



ih g* , , 9* r, i \ 

— /i B kxk'(T = ^- /i B B-er (11) 
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with fis the Bohr magneton. Thus apart from a pref- 
actor we obtain the Rashba term (Q) from Eq. (11) by 



replacing one of the k's with the electric field E. In the 
Tg vale nca b and we have two invariants for the Zeeman 
splittingcJE3 



H£ v = 2K f i B B-J + 2q f i B B-J. 



(12) 



Here, the first term is the isotropic contribution, and the 
second term is the anisotropic part. It is well-known that 
in all common semiconductors for which Eq. (^) is ap- 
plicable the dominant contribution to H% v is given by the 
first term prpopttional to k whereas the second term is 
rather smallOEil Analogous to f3\ and /?2 the isotropic 
k ■ p coupling between Tg and rg contributes to k but 
not to q. The latter stems from k ■ p coupling to more 
remote bands such, as E| and T^. 

Several authors&Eirlia used an apparently closely re- 
lated intuitive picture for the B = spin splitting which 
was based on the idea that the velocity vu = hku/m* of 
the 2D electrons is perpendicular to the electric field E z . 
In the electron's rest frame E z is Lorentz transformed 
into a magnetic field B so that the B — spin splitting 
becomes a Zeeman splitting in the electron's rest frame. 
However, this magnetic field is given by B — (v\\/c 2 )E z 
(SI units) and for typical values of E z and v\\ we have 
B ~ 2 ... 20 x 10~ 7 T which would result in a spin split- 
ting of the order of 5 x 10~ 9 . . . 5 x 10~ 5 meV. On the 
other hand, the experimentally observed spin splitting 
is of the order of 0.1 ... 10 meV. The B = spin split- 
ting requires the SO interaction caused by the atomic 
cores. In bulk semiconductors this interaction is respon- 
sible for the SO gap Ao between the valence bands Tg 
and which appears in Eq. (||). The SO interaction is 
the larger the larger the atomic number of the constitut- 
ing atoms. In Si we have Ao = 44 meV whereas in Ge 
we have Ao = 296 meV. Therefore, SIA spin splitting in 
Si quantum structures is rather small.0 

Recently, spin splitting in 2D systems has gained re-, 
newed interest because of an argument by PudalovO 
which relates the metal-insulator transition (MIT) in low- 
density 2D systems with the SIA spin splitting. Based on 
the Rashba modelo it was argued that the SIA spin split- 
ting "results in a drastic change of the internal proper- 
ties of the system even without allowing for the Coulomb 
interaction. "Q. However, as we have shown above, this ar- 
gument is applicable only to electron and LH states. The 
MIT has been observed also in pure HH systems in, e.g., 
Si/SiGe QW's.ElH As noted above, SO interact ion and 
spin splitting in these systems are rather small,Ll so that 
it appears unlikely that here the broken inversion symme- 
try of the confining potential is responsible for the MIT. 
We note that in Si 2D electron systems thep^fljective g 
factor is enhanced due to many body effects.EaS It can 
be expected that similar effects are also relevant for the 
B = spin splitting. 

In conclusion, we have analyzed the SIA spin splitting 
in 2D electron and hole systems. In 2D hole systems the 
splitting differs remarkably from its familiar counterpart 
in the conduction band. For electron states it increases 
linearly with in-plane wave vector fey whereas the spin 
splitting of heavy hole states can be of third order in k\\ . 
We have discussed consequences of this behavior in the 
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context of recent arguments on the origin of the metal- 
insulator transition observed in 2D systems. 

The author wants to thank O. Pankratov, S. J. Pa- 
padakis, and M. Shayegan for stimulating discussions and 
suggestions. 
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FIG. 1. Self-consistently calculated subband dispersion 
£±(ku) (lower right), DOS effective mass m* /mo (lower left), 
spin splitting £+(fc||) — £-{k\\) (upper right) and subband dis- 
persion £±(k\\) in the vicinity of k\\ = (upper left) for an 



\N A - N D \ = 2 x 10 1 
Fermi energy Ef- 



The dotted line indicates the 
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FIG. 2. Self-consistently calculated anisotropic subband 
dispersion £±(kii) (lower right), DOS effective mass m* /mo 
(lower left), spin splitting £f(kii) — £-(kii) (upper right) and 
subband dispersion £±(ky) in the vicinity of ku = (up- 
per left) for a [001] grown GaAs/Alo.sGao.s As heterostruc- 
ture with N s = 2 x 10 11 cm" 2 and \N A - N D \ = 2 x 10 16 
cm -2 . Different line styles correspond to different directions 
of the in-plane wave vector kii as indicated. The dotted line 
indicates the Fermi energy Ep. 
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